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Abstract 

We study the impact of the Sommerfeld enhancement on the thermal relic density of the lightest neutralino 
in the case of large co-annihilation effects with a scalar particle. The proper way of including the Sommerfeld 
effect in this case is discussed, and the appropriate formulas for a general scenario with a set of particles with 
arbitrary masses and (ofF-)diagonal interactions are provided. We implement these results to compute the relic 
density in the neutralino sfermion co-annihilation regions in the mSUGRA framework. We find non-negligible 
effects in whole sfermion co-annihilation regimes. For stau co-annihilations the correction to the relic density 
is of the order of several per cent, while for stop co-annihilations is much larger, reaching a factor of 5 in some 
regions of the parameter space. A numerical package for computing the neutralino relic density including the 
Sommerfeld effect in a general MSSM setup is made public available. 



1 Introduction 

In recent years, the existence of dark matter (DM) has 
been well established and its density measured with a few 
per cent level of precision. A recent analysis, within the 
6-parameter ACDM model, of the 7-year WMAP data [T], 
combined with the Baryon Acoustic Oscillations [1] and 
the recent redetermination of the Hubble constant 
gives g]: ^2cdm^^ = 0.1123 ± 0.0035, where D, is the ratio 
between mean density and critical density, and h is the 
Hubble constant in units of lOOkms"^ Mpc~^. 

One of the most attractive scenarios explaining the na- 
ture of dark matter is that it is composed of stable weakly 
interacting massive particles (WIMPs) , since their thermal 
relic abundance is naturally of the order of the measured 
one. However, in order to predict its precise value for a 
given model a careful calculation has to be done. In par- 
ticular, there are several effects that can alter the relic 
density of a thermal relic and which need to be taken into 
account. 

In this Letter we contribute to the study of one of these 
effects, i.e. the Sommerfeld enhancement [S]. It is a non- 
relativistic effect changing the annihilation cross section 
due to a long range force acting between slowly moving ini- 
tial particles. This effect has been studied widely recently 
(see e.g. |6j and references therein). Although its appli- 



cations were discussed for both scalar and fcrmion initial 
states, the derivation of Sommerfeld corrections for a gen- 
eral multi-state case was given explicitly only for fermions 
[7] (see also [S] for a different approaches). Here we ex- 
tend it to the case of scalar-scalar and fermion-scalar pairs 
and discuss possible applications. Finally, we present some 
numerical results for the impact of the Sommerfeld effect 
on the relic density of the neutralino in stau and stop co- 
annihilation regions in minimal supergravity (mSUGRA) 
scenario. 

2 Sommerfeld enhancement for 
scalars 

We consider an annihilation process of two particles, 
ifi and (pj , which are coupled to some light interaction bo- 
son (f), leading to a long range interaction between them. 
In the case when this interaction is diagonal (i.e. the ex- 
change of (j) does not change the particles), in the non- 
relativistic limit the spin of initial particles does not matter 
- the static force is the same for both scalars and fermions. 
This is however not true if interactions can be off-diagonal 
and intermediate particles can have different masses. Then 
due to the differences in the couplings and propagators be- 
tween scalars and fermions, the computations of the Som- 
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Figure 1: Recursive relation for the full annihilation amplitude 
including the Sommerfeld enhancement. Blob represents any 
possible annihilation process. 
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merfeld effect slightly differ. 

In deriving the effect for this case we follow the approach 
of Ref. [7] , where a general method of computing the Som- 
merfeld enhancement from the field theory diagrams was 
presented: to obtain the Sommerfeld enhancement factors 
Sij one has to solve the set of Schrodinger equations for 
the two-body wave- functions ■0^, : 



and factors A^'^ are, with </> — V, A^S indicating re- 
spectively a vector, an axial vector and a scalar: 
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and then compute 



Si- 



i'j'<p 



(1) 



(2) 



where my is the reduced mass of annihilating ipiipj pair, 
U°j contains the tree-level amplitude and the sum is over 
(possibly different) (pii (pji intermediate states and different 
interactions. Here £ = /2'mf' is the kinetic energy of 
the incoming pair (at infinity), with m'f' its reduced mass 



and p the CM three-momentum; 2Smij = rrii 



■{■ma' 



nib) is the mass splitting (for more details see Ref. [3]). 
The potential has the form: 
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where Cij^i'j' (0) are coefficients depending on the couplings 
and states involved. For incoming fermions the coefficients 
were presented in Tab. 1 of Ref. [6]. Below we give results 
for scalar-scalar and fermion-scalar pairs. 

To obtain those coefficients we use a method which is 
a straightforward generalization of the one developed in 
Ref. [7]. We write a recurrence relation for the annihi- 
lation amplitudes as visualized in Fig. [TJ Assuming that 
Smij/rrii ^ 1, in the non-relativistic limit we can trans- 
form this integral equation to the Schrodinger one, obtain- 
ing automatically the coefficients Cij^i>ji {<j)) in the poten- 
tial. 

Let's denote by a superscript S the case with two scalars 
and by F with one scalar and one fermion. In the second 
case let i and i' be fermions and j and j' scalars. Then we 
find: 
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In the limit when all the masses are equal coefficients 



Ay, 



A^ reduce to the ones which were used in Refs. 



[SKinilll]. However, in general case when the masses of in- 
termediate scalars differ, Eqs. ([SJ and ^ have to be used. 
Neglecting this fact can give rise to several percent dif- 
ference in the Cijyji{(j)) coefficients0 We stress once again 
that those results are valid in the non-relativistic limit and 
when mass splitting is much smaller than all of the masses 
involved. 

All the considerations above implicitly assumed that the 
interaction strength is sufficiently weak that the higher 
loop corrections do not alter the potential significantly. 
This is true for the weak and electromagnetic interactions, 
as well as for the Higgs exchange. However, in the case 
of strong interactions, corrections to the gluon exchange 
coming from gluon self-interactions and fermion loops may 
become important. To take this into account, following [5], 
instead of the potential ([3]) we will use one computed in 
[T^ . which in the configuration space isH 
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where ^0 = Ca — ^TpUf, n/ is the number of massless 
quarks (we choose it to be 5, since the stop co- annihilation 



^ Some of these coefficients are divergent in the limit when one of 
the masses vanish. However, in this case the non-relativistic approx- 
imation does not hold and hence these results are not valid. 

^Note that here the interaction is diagonal and one does not need 
to use coefficients II5I9I I. 
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is most important in the 0(100 GeV) region, where the 
top mass nit is non neghgible) and Euler gamma is je ~ 

0. 5772. For the case of SU(3) we have Cp = 4/3, Tp = 1/2, 
Ca = 3. For the QCD scale we take /i^ = 2mf. 

Another effect one has to take into account is the pres- 
ence of thermal corrections, as first discussed in a similar 
context in [TUj. They change the exchanged boson masses 
and in particular photon and gluon become massive, which 
introduces a Yukawa cut-off to the potential in Eq. ([TU| . 
These corrections may be also important for the mass split- 
tings, if there are nearly degenerate states present in the 
spectrum. 

There are two types of thermal effects which we need 
to include as discussed in [lOl [6] : the scaling of the Higgs 
VEV with the temperature v{T) = -T^/T^ 
where we took Tc = 200 GeV, and the contribution to 
gauge boson masses due to the screening by the thermal 
plasma; the so-called Debye mass [Tl]. For the gluon the 
screening of the plasma introduces at a leading order a 
contribution [IS] : 

ml = iNj3 + NflQ)glT^ = ^glT\ (11) 

3 Applications 

To have an idea when the Sommerfeld effect can have a 
non- negligible impact on the relic density it is useful to give 
some approximate general conditions, which have to be 
satisfied by the dark matter particle or the co-annihilating 
one: i) coupling to the boson with much lower mass ( "long 
range force" ) , ii) the coupling strength at least of the order 
of the weak coupling and iii) if the effect comes from the 
co-annihilating particle, small mass splitting between it 
and the DM. 

Those conditions are not easy to satisfy without invoking 
some new interactions. The reason is that since the DM 
has to be electromagnetically neutral, it can couple only 
to and Higgs bosons, all of which are heavy. This 

pushes up the region of possible Sommerfeld effect influ- 
ence to very large masses. Hence, without additional inter- 
actions, the only possibility is the impact of co- annihilating 
particles, which can have very different quantum numbers, 

1. e. they can even have both electromagnetic and color 
charge. In this case their annihilation cross sections can 
be altered significantly by the Sommerfeld effect coming 
form the exchange of photons and/or giuons. If such co- 
annihilating particle is degenerate with the DM, the total 
effective annihilation cross section can get large correc- 
tions. 

The co-annihilation regime is typically the only place 
where Sommerfeld effect can be significant in the minimal 
supersymmetric standard model (MSSM), or in any other 
theory beyond the Standard Model without new interac- 



tions. We will concentrate on applications to mSUGRA, in 
which the dark matter candidate is the lightest neutralino 
x5 with mass . This is because in this Letter we are in- 
terested in the effect for scalar particles and in this frame- 
work there are two extensively studied parameter regions 
involving co-annihilation with a scalar particlel^l Those are 
the stau f [19, and stop t co-annihilation regions i20i . 

The existence of the Sommerfeld effect in the stau case 
was first suggested by authors of Ref. [31] (a one-loop man- 
ifestation of this effect was also discussed in [22]), where 
they pointed out that although f+f ^ annihilation exhibits 
a strong enhancement, the f are strongly suppressed, 
and that the net result should not be very large. Here 
we will explicitly show with full numerical calculations 
that this is indeed the case, and discuss in mSUGRA its 
strong dependence on the value of tan /3. In the stop co- 
annihilation region this effect was first discussed by Fre- 
itas [9], where the QCD corrections to the bino-stop co- 
annihilations were considered, among which the Sommer- 
feld one was dominant. Below we will show our results also 
for this case0 

In order to compute the relic density, we have written a 
numerical code and implemented it into DarkSUSY [23) . 
which allows to compute the relic density with 1% accu- 
racy. 

Stau co-annihilation In this case the impact of the 
Sommerfeld effect is relatively mild. The reason is that, 
although at the freeze-out temperature the thermally av- 
eraged annihilation cross section {ucsv) is dominated by 
the contribution coming from ff annihilation, the net en- 
hancement is rather small, since as can be seen in Fig. |2] 
there are both attractive and repulsive modes present. 

Furthermore, if one enlarges the mass splitting, the co- 
annihilations become less effective and the total effect on 
the relic density gets even smaller. This has been shown in 
Fig. [21 from which one can see that the Sommerfeld effect 
introduces a several per cent correction to the relic den- 
sity, but only in the region of a very degenerate f . This 
region in most of the parameter space does not give the 
relic density compatible with the WMAP data. However, 
the impact of the Sommerfeld enhancement become more 

^In the MSSM there are other well motivated cases possible: the 
chargino co-annihilation region discussed in (6] 1161 , sneutrino DM 
scenarios [17 (however in some tension with experiment) and pos- 
sibly some small effect could be found in the "well tempered" neu- 
tralino |18| in the scenarios in which Higgs is light. Note that in 
less motivated cases, eg. when SUSY breaking scale is very high or 
when several states are degenerate with Xii the Sommerfeld effect 
gives very large corrections and is essential for reliable relic density 
calculation. 

*Our setup is slightly more general, because we do not make any 
assumption on the neutralino composition. This however do not 
introduce any significant difference in the result, since in the stop co- 
annihilation regime in mSUGRA the Wino and Higgsino component 
in the lightest neutralino is very small. 
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Figure 2: An example of the thermally averaged effective anni- 
hilation cross section, for the stau co-annihilation region. Con- 
tributions coming from different annihilation processes are in- 
dicated. The solid lines correspond to the case with Sommer- 
feld correction included, while dashed ones without it. fn the 
case the effect is too small to be visible. Parameters are: 
= 182.7 GeV, = 183.2 GeV, tan/3 = 50, ylo = and 
> 0. 

important when grows, because then the WMAP con- 
tour approaches the region where rriy^ ~ mfi the XiXi 
annihilation cross section scales as ^-nd to get the 

same (cToffw) (and the same relic density) one has to com- 
pensate with larger co-annihilation effects. This means, in 
particular, that the maximal Xi mass that can give correct 
relic density gets shifted by a Sommerfeld effect by a siz- 
able amount. In mSUGRA the value of this shift depends 
strongly on tan /3. In Fig. U we show the dependence of 
vs. the neutralino mass in the case where it is equal 
to the stau mass for its three different values. The maximal 
effect is seen for large tan /3 and drops down quite consid- 
erably when it is decreased. The reason for this is that for 
higher tan/3 the f+f^ — !■ 1^}tl^ annihilation becomes very 
efficient. On the other hand the annihilation cross sections 
of f^f^ to two leptons or quarks also grow with tan/3 but 
much slower, hence enlarging its value makes attractive 
channels dominant over repulsive ones (as in the case in 
Fig.®. 

Stop co-annihilation In mSUGRA, in the region of pa- 
rameter space where the large negative value of Aq drives 
the lighter stop to be degenerate with the neutralino, the 
results of the relic density computations are significantly 
affected by the Sommerfeld effect. 

Fig. [5] shows the change of the thermal averaged ef- 
fective annihilation cross section due to the Sommerfeld 
effect for different mass splittings between Xi and t. The 
correction is significantly larger than in the f case, be- 
cause of the strong force coming from the gluon exchange. 
When the mass splitting becomes larger two effects can be 




10-1 10-3 0.01 

Figure 3: The correction to the relic density Af2 = VLq — Qse 
coming from the Sommerfeld effect. Results for three different 
parameter points are given to show that the result is much more 
sensitive to tan j3 than m^. Note that the accuracy of the relic 
density computation in DarkSUSY is about 1%. 
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Figure 4: Effect of the Sommerfeld enhancement on the relic 
density for = ruf. Maximal neutralino mass giving relic 
density compatible with data, and its shift due to the Sommer- 
feld effect are highlighted. 
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Figure 5: An example of the thermally averaged effective an- 
nihilation cross section for the stop co-annihilation region. Pa- 
rameters are = 140.7 GeV, tan/3 = 10, Ao = -2750 GeV 
and jj, > 0. Two thresholds are visible, smeared by the thermal 
average, Xi^ s^nd it, but only the second one gets significant 
Sommerfeld correction. 
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Figure 6: Effect of the Sommerfeld enliancement on the rehc 
density for the stop co-annihilation region. For given there 
is a lower bound on due to the constraint on the lightest 
Higgs mass. 



. =255.9 GeV 
. =170.7 GeV 
= 140.8 GeV 



SI 
c? 



0.01 



0.1 



Figure 7: Sommerfeld effect on the relic density for t very de- 
generate with Xi- As before, due to the bound on the lightest 
Higgs mass the degeneracy is limited; the more the larger m^. 



seen. Firstly, the tt threshold occurs for higher tempera- 
ture, which lowers the overall impact of the co-annihilating 
particle. Secondly, also the magnitude of the correction to 
(cToffw) becomes smaller, since at higher T the typical ve- 
locities are higher and moreover, thermal corrections to 
the gluon mass are larger. 

This change in (ctcSw) can affect considerably the relic 
density of the neutralino. The results for flh^ with and 
without the Sommerfeld effect included are presented in 
Fig. inifor five different m^. One can see that the largest 
effect is obtained for parameters giving typically too small 
relic density. Nevertheless, in the region compatible with 
WMAP results, the correction can still be larger than the 
current observational uncertainty. 

The importance of the Sommerfeld effect itself is more 
clearly seen in the Fig. [71 where the ratio of relic densities 
without and with SE is presented. For very degenerate 
stops ilh^ can be suppressed by a factor of few. 

The results we presented in this subsection are in qual- 
itative agreement with those in P|. However, there are 



slight quantitative differences, for several reasons. First 
of all, in this work we were interested in the Sommerfeld 
effect and we did not compute other QCD corrections. On 
the other hand, our treatment of the Sommerfeld enhance- 
ment is more accurate, since we include not only gluon ex- 
change, but all possible interactions, and for all annihila- 
tion processes, not only for tt one. We include also thermal 
corrections which modify the masses of exchanged bosons. 

Finally, we would like to point out that the results in 
the stop co-annihilation region are subject to sizable theo- 
retical uncertainties. The reason is that since the coupling 
is relatively strong, Sommerfeld enhancement factors differ 
considerably from 1 even at high velocities. This cannot be 
however the true result, since the full quantum field the- 
ory initial state corrections in this case are not expected 
to be large. This discrepancy comes from the fact that the 
formalism used to compute the Sommerfeld enhancement 
is not valid in this regime. In our numerical calculations 
we used an approach to approximate the true corrections 
by the non-relativistic ones normalized in such a way that 
they vanish for u — 1 (a better approximation would be 
to compute the NLO vertex correction, which is however 
beyond the scope of this work)lf| To obtain more reli- 
able predictions for the intermediate regime of velocities 
0(10"^), which are very important for precise relic density 
computation, one should refine the theoretical calculations 
beyond the non-relativistic techniques used to derive the 
Sommerfeld enhancement. 

4 Conclusions 

In this Letter we studied the Sommerfeld effect for the 
scalar-scalar and fermion-scalar pairs. We gave general re- 
sults for the coefficients in the interaction potential used 
to compute the enhancement, which to the extent of our 
knowledge were not discussed in the literature before. 
Those coefficients have to be used if more than one type of 
particle is present in the computation of the Sommerfeld 
effect for a setup with scalar particles. Although, most 
applications involve only one state, there are several inter- 
esting cases in which presented results may lead to mod- 
ified phenomenological implications. Among these are: a 
general MSSM setup with heavy neutralino, > 1 TeV, 
degenerated with one (or more) sfermion, or sneutrino DM 
scenarios, again if it is degenerated with on or more slep- 
tons. It is also worth to note that new possibilities open 
up if one goes beyond the MSSM, for instance in the next 



^The normalization was done additively, i.e. the enhancement 
factors were shifted by a small constant value (always less than 1). 
Since for small velocities the enhancement is much larger than 1, 
in the non-relativistic regime this procedure does not introduce any 
significant change. This approach gives the relic density larger by 
at most 10% with respect to the case with the Sommerfeld factors 
simply extrapolated to the high velocities regime. 
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to minimal NMSSM (see eg. [24] for a review). In partic- 
ular, in the BMSSM framework [3S] one can considerably 
lower the constraints on the lightest Higgs mass, which 
would give rise to large effects coming from the light Higgs 
exchange (e.g. to Xit but even also to XiXi annihilation). 

In the second part of the paper, we presented numerical 
results for the influence of this effect on the relic density of 
neutralino in the mSUGRA stau and stop co-annihilation 
regions. For the stop case we confirm and give some im- 
provements over results of Ref. In the stau case we 
get new results saying that the effect on the relic density 
ranges from about 1 to about 7%. Typically this does not 
introduce sizable change of the which gives the correct 
relic density, however this change can be significant for the 
maximal mass, i.e. when Xi is degenerated with f, giving 
rise to about 5% shift. We observe that this effect on the 
maximal mass is largest for tan /? « 50, since then the rel- 
ative contribution to the effective cross section of f+f" 
vs. f ^ annihilation is the largest, and becomes smaller 
when tan /3 is decreased. 

For the numerical computations we have developed a 
package for DarkSUSY, which is able to compute the Som- 
merfeld effect and its impact on the neutralino relic den- 
sity for a general MSSM setup. We provide it as a public 
available tooH to be used with DarkSUSY for obtaining 
reliable relic density predictions in all possible cases where 
this effect is relevant. 
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